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I. INTRODUCTION 


The inclusion of the full pressure tensor evolution in the fluid equations of a magnetised plasma allows for a self- 
consistent description within a fluid framework of important effects that would otherwise require a kinetic treatment. 
This approach was used, e.g. to study collisionless magnetic reconnection driven by electron pressure anisotropy 
or to investigate the dispersion relation of Geodesic Acoustic Modes within the so-called Grad hydrodynamic 

equations framework recovering in this way “kinetic” effects at a much lower computational cost. In addition, a 
macroscopic point of view can help to evidentiate dynamical processes that may be difficult to identify within the 
more complex kinetic description. In particular it was recently shown that a spatial inhomogeneity in the fluid velocity 
shear can make the pressure tensor non-gyrotropic Q : in the case of a spatially localised sheared flow, the formation 
of a non-gyrotropic ion pressure tensor was shown to be mediated by the propagation in the direction orthogonal 
to the magnetic field of high frequency waves that can be interpreted as the fluid counterpart of an ion Bernstein 
branch. Thus it is of interest to examine in detail these extended fluid models in which a full pressure tensor 
dynamics is included and to evidentiate their features together with their limitations. With this aim, here we restrict 
our analysis to the case of an ion pressure tensor and investigate the small amplitude limit of the plasma description 
adopted in Ref.Q. More specifically we derive the dispersion relation of the linear waves, propagating in the plane 
perpendicular to a strong magnetic field, described by this model and compare them to the corresponding kinetic 
results obtained by an appropriate expansion of the solution of the linearised Vlasov-Maxwell system. We call these 
waves “Magneto-elastic” waves. 

The possibility of treating ion Bernstein modes within a fluid formalism, possibly extending the fluid equations to 
include high order velocity moments of the particle distribution function, may provide a convenient investigation tool 
of plasma dynamics far from thermodynamic equilibrium: for example it has been shown within a full Vlasov-Maxwell 
modelling that the excitation of a; ~ 417 ion Bernstein waves (where 17 is the ion-cyclotron frequency) in the same 
geometry configuration we consider here induces a spatially asymmetric, non-gyrotropic heating of ions @,[13. In 
a different context ion Bernstein modes are of interest in dedicated magneticall y c onfined plasma experiments on 
ion-cyclotron heating mechanisms, such as, e.g., the ALINE m and SSWIGH [l^ experiments that aim to shed 
light on the convective acceleration and heating of particles in the radio-frequency sheets, which develop both next 
to the tokamak walls because of the radio-frequency discharges and next to antennas for the ion-cyclotron heating. 
In Ref. [13 it has been experimentally demonstrated that sheared plasma flows which result from inhomogeneous 
transverse electric fields in a magnetised plasma induce the excitation of waves propagating in the ion-cyclotron 
frequency range almost perpendicularly to the magnetic field. Gonversely, the excitation of lon-Bernstein waves is 
known, both experimentally and theoretically (see Ref.fl^ and references therein), to induce the generation of 
sheared ion flows. 

Non-gyrotropic particle distribution functions are measured in the solar wind and at the solar wind-magnetosphere 
interface and correlation with the presence of shear flows has been evidenced. For example, non-gyrotropic distribution 
functions have been measured within a magnetic flux tube of ions flowing out into the upstream solar wind [T^ ■ Non- 
gyrotropic electron distributions have been observed within the electron diffusion region of reconnecting magnetic 
structures in the magnetopaus e [l8l | near X- and O-points, where the stream function of the velocity field is 
locally hyperbolic (see e.g. Ref.fl^^. 

In the following we consider the propagation of “Magneto-elastic” waves in a homogeneous Double Adiabatic 
(Ghew-Goldberger-Low, or CGL) plasma equilibrium embedded in an externally imposed magnetic field and look 
for perturbations with a wave-vector perpendicular to the magnetic field. For the sake of simplicity we take cold, 
massless electrons. We find two branches of perturbations with fluid velocities and electric fields polarised in the plane 
perpendicular to the magnetic field: a low frequency branch (LFB) corresponding to magnetosonic waves and a high 
frequency branch (HFB) corresponding to (generalized i.e., non quasi-electrostatic) ion Bernstein modes. We analyse 
their dispersion relations and compare their long wave-length behaviour to the kinetic counterparts obtained by 
linearizing the corresponding Vlasov-Maxwell system and by retaining the n = 0, ±1, ±2 ion cyclotron harmonics in 
the small Larmor-radius limit. For the sake of brevity in the following we will refer to the dispersion relations obtained 
from such a truncated kinetic Vlasov ion response as the “kinetic” dispersion relations. Qualitative and quantitative 
differences in the respective dispersion relations are evidenced. For the low frequency branch we discuss a discrepancy 
present in the literature on the derivation of the Finite Larmor Radius (FLR) corrections in the CGL-FLR fluid limit. 

This article begins with the description of the fluid model (SecHJ) and with its linear analysis in the case of ID 
perturbations propagating perpendicularly to a background magnetic field fSec lIII|) . In Sec lIVI the internal consistency 
of the model is discussed and its limitations are identified by direct comparison with the corresponding kinetic 
dispersion relations. In SeclVl] a summary of the results obtained is presented and the possibility of describing 
higher ion cyclotron harmonics (|n| > I) within a fluid framework is noted. Detailed calculations are presented in the 
Appendices: in Appendix|^a brief derivation of the two-fluid model equations is given, in Appendix [B] the limit of the 
Double Adiabatic closure is described and in Appendix [C] the derivation of the FLR limit of the magnetosonic branch 
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is detailed. We recall in this context that the Double Adiabatic closure (see e.g. M) describes the low frequency, long 
wavelength dynamics of a magnetized plasma with different temperatures, and obeying different equations of state, 
in the direction of the magnetic field and in the plane perpendicular to it. It assumes temperature isotropy in this 
plane. This closure applies to plasma conditions where heating and cooling processes and/or approximate integrals 
of motion due to the presence of particle adiabatic invariants affect the parallel and the perpendicular directions 
differently while it assumes that the smallness of the particle gyroradii and their fast gyration are sufficient to allow 
for a fluid isotropic description in the perpendicular plane. For the perpendicular propagation considered here only 
fast magnetosonic waves could propagate if the Double Adiabatic description were applied to the perturbations of the 
equilibrium. 


II. FLUID MODEL WITH FULL PRESSURE TENSOR EVOLUTION 


We consider a single-fluid MHD model in which the contributions of the electron pressure and of the electron inertia 
are neglected and the ion pressure evolves according to the full pressure tensor equation. This model can be obtained 
from the two-fluid, full pressure tensor equations derived in Appendix by adding in the standard way the ion 
and electron momentum equations fEa. (IA3l) l and by using the quasi-neutrality assumption n® ~ n® for a hydrogen 
plasma. Thus the magnetic force and the spatial derivatives of the ion pressure tensor are the only forces acting on 
the plasma (Eq.^), while the electron momentum equation appears in the single-fluid system in the form of the 
Hall-MHD Ohm’s law in the induction equation (Eq.(IS|)). Furthermore we disregard the ion heat flux, diQijk = 0 in 
the ion pressure tensor equation (Ea. (IA6l) l. In the present analysis this simplifying assumption is made as a closure 
condition whose consistency with the truncated Vlasov model in the low-/3 limit will be shown later fSec lIVI) . Here 
P is the ratio (suitably defined, see Eo. lflTl) ! between the ion pressure and the magnetic pressure. We note however 
that, as stressed in Ref. (2ll| . the comparison with the kinetic dispersion relation may be made exact in the context of 
the low frequency limit of magneto-sonic waves also at /I ~ 1 if the heat-flux contribution is consistently retained (see 
SecIVland Appendix ICl) . 

Introducing the ion cyclotron frequency D and the unit vector components bi = B^/B, the single-fluid equations read, 
with standard notation, 

dtn -I- V • (nu) = 0, (1) 


dtU + u ■ 'Vu = n — X b — 
ne 


V n 

mn ’ 


( 2 ) 


dtU + V ■ (uU) + n ■ Vu + {n ■ Vu)^ 


( 3 ) 


-D(n X b-f (n X bf) = 0. 


The apex “T” indicates matrix transpose. Note in passing that Eq.([3]) ensures that the pressure tensor remains 
positive definite over time if it is positive definite at t = 0. 

In the adopted model the displacement current is neglected in Ampere’s law (Eas. dATI) '). as consistent with the quasi¬ 
neutrality assumption and the restriction to phase velocities smaller that the speed of light thus leading to the MHD 
definition of the current density, 

J = £vxB. (4) 

Then the induction equation coupled with Ohm’s law is 

u -^ X sj’ . (5) 

The third and the fourth terms in Eq.Q arise from the action on the pressure tensor 11 of the spatial derivatives of the 
velocity field u which we write in tensor notation as the sum of a strain (that is, of compression and of compressionless 
shear) and of a vorticity term as 



duj 

dxj 


1 

3 



compression rate 



/ dui duj \ 2 / duk \ 
\ dxj dxi J 3 \ dxk ) 

compressionless rate of shear 



1 / dui duj'\ 

2 yidxj dxi) 

fluid vorticity 


( 6 ) 


strain tensor 









4 


III. LINEAR ANALYSIS OF THE FLUID MODEL 


We linearise Eqs.dTJSl) around a double adiabatic equilibrium II?- = + {P^ — P^)bibj (see ApplB]) with uniform 

density no, no velocity = 0), uniform magnetic field = BqBz. A double adiabatic equilibrium is the most 
general, spatially uniform configuration compatible with the pressure tensor equation (Eas. (IBlllB2l) l in the absence 
of a shear flo w 1^ . When the latter is included, a class of spatially inhomogeneous, non-gyrotropic equilibria can 
be found l24||. Here we consider perturbations with a wave-vector perpendicular to the equilibrium magnetic 
held, k = kex, since this makes it possible to decouple the dispersion relation of modes with electric and velocity 
held components perpendicular to the equilibrium magnetic field BqBz from the modes with parallel components. 
This simplifies the study of the propagation of perturbations that make the pressure non-gyrotropic in the plane 
perpendicular to BqB^- Note however that the curl of the Hall term in Eq ([5]) vanishes in this geometry. 

Labelling perturbed quantities by a ~ we set 




It is convenient to introduce the characteristic Alfven velocity (ca) and the ion thermal velocity (vth) as 








47rnom “" m nom 

We then introduce the ion skin depth (d^) and the thermal ion Larmor radius (pi). 


( 7 ) 


( 8 ) 


df = 


In this model in which only the ion temperature contributes to the plasma pressure, the scale separation between di 
and Pi is measured by the magnitude of the P parameter, 


'^th 


(9) 


cl dr 


Linearising Eq.([2]) in this geometry we find the equations for the in-plane velocity components. 


' k\ 2 - , k Uxx 

Ux — I I C^Ux -j- , 

w y w nom 


k H 


u„ = — - 


xy 


uj nom 

which are coupled to the equations for the perpendicular pressure components, as obtained from Eqs.Q, 


w, 

n 


H — IP' 


0 kux 


n 


2in. 


xy, 


—H — pO I j(f\ _ fr ') 

— p_l ^ w l-yxiyy l-^xxJ, 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


U) ~ n kUx 

_TT — pO ^ o„TT 

^ yy ~ -L ^ 


(15) 


and to the equation for the parallel pressure component. The latter evolves because of the plasma compressibility, 

wH^z = P°^kux. (16) 

After a few algebraic steps the linearised system equations can be cast in the form [M] ■ u = 0, with the dispersion 
matrix given by: 


/ krcl + v^h) I _ k 


[M]^ 


2 2 
Cth 


,11 krr 


' LU (402 _ 


H k^vfh \ 

2(402 _ ^2^1 ^ (402 _ ^2^ 

1-k' 


^th 


2(402 -a;2) 


(17) 


V 
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A further set of closed equations is provided by the parallel component of the velocity, Uz, which is decoupled from 
the perturbations (■Ua,,'Uy,0) but is coupled to the linear equations for Il^z and Tlyz, 


k IIt 


Uz = -- 


u! nom 


(18) 


^f\ - 


*n„ 


(19) 


‘^liyz = -ihxz- ( 20 ) 

This set provides the dispersion relation of a spurious mode which is unphysical having no counterpart in the Vlasov 
description; 


1 + 




2{n^-uj^) 


= 0 , 


( 21 ) 


(the differences with respect to the corresponding branch obtained from a kinetic treatment are evidenced in Sec IlVI) . 
In fact this mode would correspond to a hydrodynamic mode with no perturbed electromagnetic fields and thus cannot 
be realized in a collisionless plasma. This feature originates from the simplified description of the parallel dynamics 
adopted in the present fluid model which, for k T Bq and in the limit of cold mass-less electrons, does not make 
it possible to balance a parallel electric field component but it does not affect the modes with perturbed velocity 
perpendicular to the equilibrium magnetic field treated in the rest of this paper. We call the modes described by 
Eqs. (umii) “ magneto-elastic modes”, as they arise because of the pressure anisotropy which couples the electromag¬ 
netic perturbations described by Ea. dTTl) to elastic-type perturbations related to the Uy component lEa. dT^ l through 
the off-diagonal component Bxy. The evolution of IIz^, though not directly contributing to the dynamics of u^_, is 
coupled to it (see Ea. dTBl) ') because of the compression related to the velocity component Ux fSec lIII AT) . 

When comparing the full pressure tensor dispersion relation with the isotropic MHD model where the magnetosonic 
branch has dispersion relation = k'^{c\ -I- c^), with Cg sound speed of the plasma, we need to recall that, since 
the electrons are taken to be cold, here the sound velocity is replaced by Vth of Eq.®. In particular, the factor 2 in 
the definition of corresponds to a polytropic index Tx = 2, as follows from the double adiabatic equation for P± 
lEa lBTI) . when Bq - k = 0. 

It is convenient to write the dispersion relation of the modes u = u± as 


UJ 

112 


-4- 




112 


- k^ 


- 2k^P^ 


= 0 , 


a;2(4112 — uj) 


( 22 ) 


which describes the two branches 



= 2 + k^ 




2 ’ 


(23) 


where, consistent with the notation of Ref.Q, the higher frequency branch (HFB), corresponding to the root with 
the plus sign in front of the square root, is denoted by Wh while the lower frequency branch (LFB), corresponding to 
the minus sign, is denoted by uii. The behaviour of these dispersion relations is exemplified in Fig[Tl left frame. For 
UJ < 20, only the LFB propagates. It corresponds to a fast magnetosonic wave, corrected by the presence of the full 
pressure tensor evolution. From a local expansion of the matrix [M] we see that this branch is compressive and has 
a purely transverse electric field at low frequencies (w/H <C 1), with eigenmodes of components 


{Ux,Uy) 



n k'^Pi\ 


{Ex, Ey) — 


^ /I ’ ^ 

uji 4 


(24) 


while it becomes approximately right-hand circularly polarised, {ux,Uy) ^ i{Ex,Ey)c:i (IjO ^-t the crossing of the 
“resonance” uj ~ 211. At a; = 211 the HFB appears at fc = 0: its dependence on k is initially flat and grows linearly for 
larger values of k. This HFB is related to the n = 2 “generalised” (i.e. not quasi-electrostatic [1^, [Ig) ion-Bernstein 
mode which is found in a kinetic description (see Sec lIV AT) . For kdi <C 1 this mode is left-hand circularly polarised 
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FIG. 1. Dispersion relations (left frame) and group velocities versus kdi (right frame) of the LFB {coi and vi, dashed lines) and 
of the HFB {toh and Vh, solid lines) obtained from the fluid model for different values of vth/cA = ps/di. 


with {ux,Uy) ~ i{Ex,Ey) ~ (1, —i), and maintains approximately the same polarisation for kdi ~ 1 (as deduced by a 
local expansion of the eigenmode equations around the resonance). This branch is not present in a double adiabatic 
description and it arises within a fluid description because the full pressure tensor equations allow for non-gyrotropic 
pressure perturbations in the x-y plane. The lower frequency bound of this branch, lj = 217, follows from the fact that, 
because of the equilibrium magnetic field, these components rotate in the x-y plane at twice the cyclotron frequency. 
The group velocity of the two branches, normalised to the Alfven velocity, can be expressed as 


Vh,l _ 1 I I ^ ] 

Ca uJh,i 1 2 df B j 


(25) 


where the ± sign in front of last term corresponds to the high and to the low branches respectively and 


A = -l + k^ 



4dj ’ 


(26) 


B + (27) 

Its behaviour is exemplified in Figll] right frame, for some values of di and pi. 

Note that in the limit pi —>■ 0 (cold ion limit, vanishing pressure tensor) the two branches come close to crossing each 
other (see left frame of Figll]) and that for pi = 0 the HFB disappears (in formal terms the solution txP '— 4 = 0 
of Eg. (1^51) is cancelled by the corresponding resonant term in the denominator of Ea. (0^ ). 

Ordering pi ^ di, in the long wavelength limit {pf ^ df, /c —>• 0) we obtain 

uf ^ k\cl + vl„), u^l 40^ + (28) 

whereas in the opposite case of very short wavelengths (p| ~ df, fc —>■ oo where, however, a fluid-type description is 
not expected to hold as mentioned at the end of this section) we have 


OJ 


2 

I 


^th 


-b -V, 


'^th 


(29) 


Note that in this formal limit the off-diagonal terms of the matrix tend to zero and thus, as expected in the case 
of short wavelengths, the two branches separate into a purely electrostatic (w^ ~ + (3/2)U(,^]) and a purely 

electromagnetic {ujf ~ fc^U(^/2) mode. 

Eor comparison with the CGL-ELR model discussed in Sec lIII Bl it is interesting to consider the first order corrections 
to the LFB, obtained in the small Larmor radius limit by performing the ordering <C and fc^pf <C |4 —fc^d^j. 
We find 


c 


2 

a 



\ 

2{4-k^d^)) 


^th 


1 - 


4 - 2fc^df \ 

(4-fc2d|)V 


2 7 2 
uji ~k 


(30) 
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The inequality k'^pf <C |4 — k'^df\ follows from the fact that a small expansion is not valid in the vicinity of the 
near crossing between the two branches mentioned above. 

We conclude this Section by noting that, when compared to the corresponding dispersion relations computed from 
the truncated Vlasov expansion in Sec lIVI the range of validity of the dispersion relations represented in FiglHwill be 
restricted to values kpi < 1. Since kpi = kdi\/P^ this condition becomes more restrictive with increasing values of /3 
when expressed in terms of di. Beyond the interval kdiy/P < 1 it is not possible to provide an interpretation of the 
fluid hierarchy closure in terms of an expansion of the Bessel functions (see Sec lIVIfVl) . Clearly the fluid description is 
expected to fail for kdi ^ 1, with its range of validity depending on the branch considered. 


A. Compressional effects and fluctuations of 

Both branches induce fluctuations in and II^^. These fluctuations are coupled through the relationship H^z/B^ = 
P|*|/i3o, since linearisation gives 


,01', 

po “ w ’ Bo~ uj ' ^ ’ 

This polarisation is coherent with an isothermal closure for the parallel “temperature”, Ilzz/inom)^ where n/uQ = 
kux/^ is obtained from linearisation of Eq. o, since for k = kxS^ the double adiabatic equation for P|| (Ea. (IB6l) l 
corresponds to a polytropic with index Tn = 1. 


B. Dispersion relation in the limit of a CGL closure with first order FLR 


For further comparison we consider the well known CGL fluid model with first order FLR corrections [27 h 32|| 
recently re-derived in Ref. |23j by starting from a full-pressure tensor model and by assuming uj/fl ^ kvth/^ ~ e 1. 
By construction this model projects the linear system of Eq. (ED onto the LEB, since it relies on the small frequency 
approximation w/il 1 at the basis of the CGL closure. However, as first noticed in Ref.fs^ and in Ref. [2^, it 
fails to describe this magnetosonic branch correctly for the purely perpendicular propagation that we are considering 
here. This fact was further pointed out in Ref.[2l| while discussing how the CGL model with first order FLR 
corrections may lead to an incorrect description of transversely propagating, low-frequency magnetosonic solitons. 
Here we remark that the failure of the fluid description of fast magnetosonic waves in the CGL-small FLR limit may 
be understood by noticing that for ■ k = 0 the CGL closure with first order FLR corrections does not take account 
of the correct polarisation of the magnetosonic waves which requires Uyjux ~ kdi “C 1. For • fc = 0, indeed, the 
maximal ordering Ux/vth ^ Uy/vth ~ , which the CGL-first order FLR model relies on, is inconsistent with the 


neglect of corrections in Ea. dHl) because from Ea. dT^ it follows that Uyjvth 
wave-length limit, ^ <C 1, of Eg. ll^Op . 


e. On the other hand, the long 


2 ; 2 
uji — k 


k'^pf \ 2 

1-^1 


1 - 


16 


(32) 


Ca 


kCg 


1 - 




Pi. 


1 - 


Ppt 


(33) 


is directly obtained from the set of Eas. dlHIlS]) after ordering w/H ~ kvth/^ ^ Uy/vth ~ and Uxivth ~ The 
relevant calculations are detailed in Appendix o while in Sec II V A] there is the dispersion relation lEa. (l53ll l to be 
compared with the small-/3 limit of Eg. (15^ . Here we recall that by linearising the GGL-first order ELR set of 
equations for perpendicular propagation we would obtain instead 


2 ; 2 
uJi =k 


[ 1 + 


2 „ 2 \ -| 


16 


(34) 


with group velocity 


n 

Cn 


kCg 

UJl 


1 + 


d? 


k^pt 


(35) 


















4 
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1 
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kdj 
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FIG. 2. Dispersion relations (left frame) and group velocities versus kdi (right frame) for the LFB as obtained by retaining 
small FLR corrections to the CGL closure for different values of Vthjca = pijdi = the frequency w; from Ea. dMll and 
the group velocity vi from Eg. (1351) are represented by dashed lines. For comparison with the full pressure tensor model, the 
curves of toi and vi obtained from Eas. (l23l25ll are also represented (solid lines) for the same parameters. For the red dashed 
curves {fd = 10“^) the condition kpi < 1, required by the small-FLR expansion, is fulfilled over the whole interval. For the 
black (/?=!) and blue (/3 = 1/4) dashed curves the formal ranges of validity kdi < 1 and kdi < 2 are delimited by the black 
and blue dotted vertical lines, respectively. Even if the above intervals also dehne the range in which curves with equal colors 
should be compared, no formal restriction on the wavelength range is assumed for the solid lines which are computed from the 
full pressure tensor fluid model without expanding for kpi < 1 . 


This latter dispersion relation, first given in Ref. and obtained by neglecting the leading order perpendicular 
(to B) heat-flux contribution as for Eos. (13211551) . misses the FLR correction to the Alfvmiic contribution of Eg. (15^ 
because of the incorrect polarisation (Appendix ICl). In addition the FLR correction to the thermal velocity appears 
with the opposite sign. While the term should be neglected in the small-/? limit k^Pi k^d^, in which 

Eq. (15^ coincides with the kinetic counterpart obtained by consistently neglecting the heat-flux gradient contribution 
(Sec II V AT) , in Eg. (1551) it represents the only dispersive effect, which is of the same order as the neglected perpendicular 
heat-flux term. 

In Figs.(I2]) the dispersion relation for the LFB (left frame) and the dependence of the group velocity on k (right 
frame) at relatively small values of kpi are compared for the solution obtained in the fluid model with a full pressure 
tensor equation (Eas. (l23l25L solid lines) and in the CGL-FLR model (Eos. (l34ll551) . dashed lines) for different values 
of pi/di- The range of values of kdi in Figs. ([2]), corresponding to values of pi that extend beyond the limit kpi ^ 1, has 
been chosen in order to highlight the differences between the full pressure tensor model and the CGL-FLR limit. This 
applies in particular to the cases pi = O.Sdi (blue curves) and pt = di (black curves) in which the small FLR limit is 
often considered, especially for astrophysical applications (see e.g. the discussion in Ref. [2^). Even if the differences 
may appear negligible in the wave-length range kdi ^ 2 and kpi -C 1, a very different behaviour is displayed for 
kdi ^ 1, even if k'^pf —>■ 0 (red curves). These differences become remarkable for the group velocities, as the small 
FLR limit obtained from Eq. (1501) gives a group velocity, see Eq. (1551) , that decreases with increasing values of kdi , 
a behaviour in agreement with the kinetic result [2l|, [2^ (cf. Eg. (155)) next) and in contrast with the increase 
described by Eo. (l35l) . 


C. Eigenmodes in the long wave-length limit 

In the long wavelength limit kdi ^ 1 the eigenvectors V’'{k,uJi) and V^{k,ujh) in the perturbation space {ux 
BxyVth/P±, {Bxx — ^yy)'>Jth/P±} of the LFB and of the HFB are respectively 


W= ul<l 


2^2 


n k^pf n k^p 




4 ’ kvth 4 


2„,2 

th 


k^ / 2iof-k^ v. 
Auji I kvth 


(LFB), 


(36) 


yh. ^ 


I5 b i'^th , 5 ^ 


uJh 1 f2ujl-k^v^f^ 


fc ’ n y kvth 


(HFB). 


(37) 
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where the relative sign between the pressure perturbations and the velocity perturbations depends on the propagation 
direction i.e., on the sign of the phase velocity. Note that, because of the presence of the two branches and of their 
different propagation properties, an initial incompressible perturbation of the form V = {0, uo{x), 0, 0} can generate 
compressible flows since the difference in the propagation velocities removes the cancellation that makes Ux {x) = 0 at 
the initial time. This feature is central to the mechanism of shear-induced agyrotropy discussed in Ref.@ and cannot 
occur in the CGL-FLR limit where only the uji mode is present. 


IV. KINETIC DESCRIPTION OF THE LOW AND OF THE HIGH FREQUENCY BRANCHES 

While the low frequency magnetosonic branch has its counterpart in the MHD plasma description, this is not the 
case for the high frequency branch which cannot be described by the MHD equations. In order to characterise the 
origin and nature of this latter propagation branch more precisely, in this section we consider for both branches the 
corresponding kinetic dispersion relations as obtained from the solution of the linearised Vlasov-Maxwell equations. 
We assume again purely perpendicular propagation. 

First we note that in order to account for the formation of a non-gyrotropic pressure tensor in the x-y plane (i.e., in 
the plane perpendicular to the uniform magnetic field = EqE^), in the solution of the linearised Vlasov equation 
we need to include the m = ±2 terms in the expansion of the perturbed ion distribution function into harmonics of 
the gyration angle a in phase space. The combination of the contributions of these terms to the kinetic dispersion 
relation leads to the — 417^ resonances that we have encountered in the fluid treatment. For a similar reason, 
having neglected in the fluid analysis the ion heat flux diQijk, we need not include the m = ±3 (and higher order) 
terms in this kinetic treatment. 

Thus we compute the ion contribution to the permittivity tensor Kij by retaining the n = 0, n = ±1 and n = ±2 
terms in the summation over the cyclotron harmonics (see e.g. Ref.js^, p.404-405). In addition, in order to make the 
comparison with the fluid treatment meaningful, we take the small Larmor-radius limit by retaining contributions up 
to the power ^ k'^dfpf while neglecting those in ~ This means that in the elements of the permittivity tensor 

we have kept up to the linear contributions in {kpi)^/2 and that the comparison with the fluid derivation is restricted 
to values (fc^pf)/(/c^df) = /3 <C 1. Clearly this expansion procedure does not allow us to recover in the low frequency 
limit the FLR correction to the k^v^f^ in Eg. (15^ . On the other hand, it elucidates the inconsistency of the dispersive 
effects in the CGL-FLR dispersion relation Ea. dMl) . These specific points will be considered later in this Section. 

We compute the electron contribution to the permittivity tensor taking the limits —>■ 0 and —>• 0, with 

Qg = eBQ/{meC) and pe the electron cyclotron frequency and thermal Larmor radius respectively, as consistent with 
the zero mass (cold) limit assumed in the fluid treatment. Thus we retain only the off-diagonal contributions to the 
permittivity tensor of the E x JB-drift which are given by the low frequency limit (w^/Qg <C 1) of the n = ±1 terms in 
the summation over the electron cyclotron harmonics. By adding the ion and the electron contributions in the limits 
described above we obtain 


i^ii = 1 - oji 


l-k^pl/2 2 


k^p^ 


LU 




2(a;2-4Q2) 


2\ ’ 


n=±l 


n=±2 


(38) 


Kt2 = —/C 21 = —i 





n—dzl 


n=i2 


(39) 


2 1 ^k'^p‘f/2 2 

A 22 = I — -^-^- — - U}„ 


2„2 


k^p: 


Q2 


22 - Q? 


n—0 


n=±l 


2(a;2-4Qf)^ 


n=±2 


(40) 




A33 = 1- 


2Q2 


Pp2 




2(a;2-Q2)- 


n=±l 


(41) 


n—0 
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where Wpg/fie = The order n of the cyclotron harmonic from which each term comes is labelled in the 

underbrace. Note that the cancellation between the electron and the ion E x B contributions to the plasma current 
leads in the low frequency limit to a partial cancellation between the first two terms in if 12 , and if 2 i- 

Neglecting the displacement current, as consistent with the fluid model used above, i.e., by taking <C c^, we 
write the normalised dielectric tensor Dij = ( 0 ^/a;T){ifij — — kikj)c'^/uj^} as 


/ - k^v\j2 


\D] = 


+ 


^th 


n 


n / (jj 

—i— 


2(4122 _a;2) 
2 




^th 


<^th 


+ 


^th 


V 


w V 12^ — a;2 4122 _ ^2 

0 

The vanishing of the coefficient 

D 33 = 


uj \ 122 — a;2 4j^2 _ ^2 


OJ^ 


^th 


122 _ ^2 
0 

n^ + k^{cl-vfj2) 


2(122 -a;2) 

yields the dispersion relation of the so-called Ordinary mode, 

9 9 ' 

uj^ = n^ - 


'^th 


2(1 +Pd?)’ 


^th 


0 


\ 


2(4122 -P) 


0 

D 33 J 


(42) 


(43) 


(44) 


which corresponds to electromagnetic perturbations with E parallel to Bq. This mode has no counterpart in the fluid 
description in Sec. IIIII where, as already mentioned, the “hydrodynamic” root related to Uz in Eq. (II3 does not apply 
to a collisionless plasma. 

The dispersion relation of the transversely polarised modes, obtained from the vanishing of the determinant of 
Ea. (H^ by consistently neglecting the ~ k'^pf contributions while keeping the ~ k'^d^p} terms, is 


4- 


UJ 

IP 


w 

122 


- k^df (1 - 


fc 2 p 2 


fc 2 p 2 


fc^d? ( 1 — 


UJ 

122 


T = 0. 


a;2(122 — a;2)(4122 — uj"^) 

Its roots, whose behaviour versus kdi is sketched in the left frame of FiglHl are (same notation as in Ea. (|23)) ') 


(wT = 2 +^ ± \/^ - +^ (<«+ 


The group velocities, represented in Fig. ([3]), right frame, are given by 

Vh,i _ Ji ± 

Co ujh,i 1 2 b] ' 

where 


A = -1-2 


d? 


1 + 6 


d? 


(45) 


(46) 


(47) 


(48) 


B = 


h2p2 

4 - 2 Pd? - 4 fc 2 p 2 + (fc 2^2 


6 Pp?). 


(49) 


As in the fluid description, see Eos. (12211271) . in the limit pi = 0 the HFB vanishes as indicated by the cancellation of 
the 412^ — uj^ = 0 terms in the numerator and in the denominator of Eq. (|45|) . 

The LFB and BFB described by the roots of Ea. (H51) are respectively the magneto-acoustic mode and a generalized 
ion-Bernstein mode. In the quasi-neutral limit considered here this mode is not quasi-electrostatic, as consistent with 
the general dispersion relation discussed in Refs. [ 2 ^, (^ . 
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FIG. 3. Dispersion relations (left frame) and group velocities versus kdi (right frame) of the LFB (cj; and vi, dashed 
lines) and of the HFB {ujh and Vh, solid lines) obtained from the truncated Vlasov-Maxwell system for different values of 
Vth/Ca = pi/di= y/P < 1. 




FIG. 4. Dispersion relations (left frame) and group velocities versus kdi (right frame) of the HFB uJh obtained from the 
truncated Vlasov-Maxwell system (red lines) and from the fluid model (black lines). Both the HFB (solid lines) and the LFB 
(dashed lines) are represented. The value Vthjca = Pi/di= y/p = 0.01 has been chosen as representative for the <C 1 limit 
in the interval of wave-numbers considered. 


A. Comparison of the fluid and kinetic dispersion relations 


Because of the assumptions made in the derivation of the kinetic dispersion relations (Eas. lHSlH^ ). a meaningful 
comparison with the fluid counterparts (Eas. (l22ll27)l ) is possible only for small values of i.e., for relatively small 
values of the ratio pi/di= y/p, and for a restricted range of values of kdi. This comparison is shown in Fig.(|3]). 


First we note that in this small limit the HFB described by Ea. (H51) has, near the resonance to = 2il, the 
factorized expression 


Pi (fc^d2)(3fcHV2-4) 

H ^ dy^ df 2(4-Pdf) 


We see that the dependence of dw = w — 2H on pf appears only through a multiplicative factor whereas the values 
of kdi corresponding to its minimum (zero of the group velocity) and to the second crossing of the resonance at 
fed, = 2v^, are independent of pijdi. On the contrary in the fluid model the HFB remains above the uj = 2Q 
resonance for all k 0. This different behaviour is made evident by the opposite spatial dispersion of the HFB for 
—>■ 0 in the two models where we find 


Vh^± 


kvjh 

2n ’ 


col c 4H2 ± k\l, 


(51) 
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where the plus sign refers to the fluid model and the minus sign to the kinetic treatment. This opposite behaviour is 
also shown by the first order correction to the ratio of the perturbed pressure tensor components, 

n,,/n,,~-(i±fcV?), (52) 


where again the plus sign refers to the fluid model. The perturbed pressure components are computed in the fluid 
model directly from the linearised pressure tensor equations and in the kinetic treatment from the second order 
velocity moment of the perturbed distribution function (Ref.js^, Eq.8.10.8). 

Secondly we note that while in the fluid model the LFB magnetosonic branch crosses the uj = 20 resonance at 
a value which in the limit <C 1 is kdi ~ 2-\/2 (as obtained from Ea. (l23)l 'l. in the kinetic treatment this branch 
stays below this resonance. More importantly, having retained only the corrections in both the diagonal and 
off-diagonal matrix elements of the dielectric tensor, the low frequency limit of the kinetic dispersion relation (Eg. (1461) 1 
can account consistently only for the ELR corrections to in the dispersion relation of the magnetosonic branch. 

The magnetosonic dispersion relation obtained by disregarding the ^ k^pf terms in the expansion of the LFB of 
Eq.dlSl) is 


Wz 


k^ 



+ 


(53) 


This equation agrees with the small-/3 limit of the kinetic result obtained in Refs. [2l|, [2^ and coincides with the 
consistent limit of Ea. (|32ll . 


V. ROLE OF THE HEAT-FLUX IN THE LONG WAVE-LENGTH LIMIT OF THE LFB FOR /3 ~ 1 


Neglecting ^ k'^pf contributions in the power expansion of Ea. (H51) allows us to obtain the same LFB dispersion 
relation (l53)) from the truncated Vlasov-Maxwell system and from the fluid model, closed so as to include the full 
pressure tensor equation without the divergence of the heat flux. This however does not imply that the heat-flux 
can always be consistently neglected even when only terms that scale with the perpendicular wave number as k^ are 
included. 

For j3 of order unity and for purely perpendicular propagation, terms of order must be retained in 

the dispersion relation of the LFB in the long wave-length limit. In the fluid description a term of this order arises 
from the heat flux in the limit u! and needs to be retained for a consistent description of the LFB, see Appendix 
o This feature of the LFB was first pointed out and extensively discussed in Ref.[^. In the kinetic derivation a 
term of the same order arises from the n = 0 contribution to the K 22 element of the permittivity tensor which, for 
uj/rt <C 1 , requires us to retain an additional ~ term in the expansion of the Bessel function since this term 
contributes to the dispersion relation of the LFB to the same order of accuracy, in powers of (w/fl)^ ^ k^Pi, as the 
other contributions to Ea. (l40l) . In doing so Ea. (l40l) becomes 


K 22 = 1 - 


^th 


II 2 


1 - 


’pi 2 1 3fc^p?/2 2 

J ~ a;2 -02 


2 


k^p: 


2(cc2 - 40f) 


n—0 


n=±l 


and the D 22 element of matrix (1421) is rewritten as 




D 22 = — - 


oi + vtAi--k^pl 


02 - ik^vU^ 


n=±2 


^th 


w2 02 -a;2 2(402 -^2)’ 

This leads to the correct extension of Eg. (1551) including ~ k'^vf/^k'^pj terms, which reads [2l|, [2 


2 7 2 

uJi :^k 


k'^p^ 


1 - 




(54) 


(55) 


(56) 


Finally we note that the |ri| = 1 harmonics are absent from the dispersion relation of both the fluid and the 
“truncated” Vlasov analysis. In order to recover the |n| = 1 harmonics in the dispersion relation obtained from the 
truncated Vlasov derivation terms proportional to (fc_Lpi)^ must be retained and the limit to w —^ 17 and (fc_LPi)^ —>■ 0 
must be performed in an appropriate order. Correspondingly, it can be shown that the |n| = 1 harmonics are recovered 
within the fluid description if the divergence of the heat flux tensor is retained in the equation for the time evolution of 
the pressure tensor. Consistently with the truncated Vlasov formulation this contribution is proportional to (k±pi)'^. 
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VI. CONCLUSION 


The aim of this article is to characterise the linear modes that can propagate in a homogeneous plasma in the 
direction perpendicular to an externally imposed magnetic field within a two-fluid description that uses a simplified 
model for the cold electron fluid and a full pressure tensor treatment of the ion fluid. The perturbed velocity fields 
of these modes lie in the plane perpendicular to the magnetic field. This model description in its nonlinear version 
has been used used in the investigation of the generation of a non gyrotropic pressure tensor Q . It is thus important 
to have an a priori understanding of the physical properties of this model and of its limitations, in particular since 
it has been examined in full detail only in the context of its equilibrium configurations [^ . A first step in this line 
is to determine which kind of waves this model describes and in particular how well they correspond to the results 
obtained by solving the Vlasov equation. Such a comparison is more easily performed under the homogeneous plasma 
conditions adopted in the previous sections. 

More specifically we have compared the linear dispersion relations derived from this fluid model with those obtained 
from the solution of the Vlasov equation in a magnetised plasma by retaining, in the appropriate expansion in the 
thermal ion Larmor radius, the contribution of the n = 0, ±1, ±2 ion cyclotron harmonics. The latter harmonics 
are required in order to account for the dynamics of an anisotropic, and in general non gyrotropic, pressure tensor. 
This comparison shows that plasma kinetic dynamics involving the second cyclotron resonance can be mimicked, in 
the appropriate small Larmor radius limit, by a fluid description provided a full tensor pressure dynamics is retained. 
Even if the agreement between the two descriptions is not exact in some cases, as explicitly mentioned in our analysis, 
the overall framework is coherent enough to make the extended fluid model viable and useful when spatially extended 
inhomogeneous configurations are considered. In such a case fully kinetic simulations are exceedingly expensive in 
computational terms. 

Besides, this fluid approach makes it possible to identify, more easily than with a full kinetic description, some features 
of the normal modes propagating in the system, such as the way the dynamical plasma response to perturbations 
determines the polarisation of these normal modes. 

The inclusion of ion thermal fluxes in the plane perpendicular to the magnetic field would extend this description 
making it possible to recover the n = ±I cyclotron harmonics that do not appear in the dispersion relations derived 
above and to include in principle higher harmonics together with higher order terms in Bessel function expansion. 
In addition, in Sec. El we have shown that, while corrections arising from higher order terms in Bessel expansion of 
the n = 0 cyclotron harmonic appear to be negligible at high frequencies w > 17 for small k, they become important 
when considering FLR-corrections to the magnetosonic dispersion relation at non-negligible values of /3 as evidenced 
in in terms of the contribution of the perpendicular heat flux. In this context of comparison between fluid and 

kinetic frameworks, we recall that the dispersion relations of magnetised plasma modes such as kinetic magnetosonic 
modes and, using a different approach, Bernstein waves, have been reconstructed from a fluid-like analysis in Refs. 


More specifically we have shown that the fluid model describes two wave branches that involve velocity perturbations 
perpendicular to the magnetic field: a low frequency (LFB) and a high frequency (HFB) branch with dispersion 
relations ojf ~ k'^ic^ + and ~ AD? + respectively in the long wavelength limit kdi 0. While the 

LFB corresponds to the fast magnetosonic waves propagating in a cold-electron plasma, where the thermal ion speed 
replaces the usual sound-velocity, the HFB finds no fluid counterpart in the isotropic pressure limit. The origin of 
the HFB can be understood by comparison with the kinetic dispersion relation in the small FLR limit ISec lIVI) . from 
which we recognise that the HFB corresponds to a generalization of the 217 ion-Bernstein mode. 

Finally we underline that the so called CGL-FLR fluid limit, obtained by expanding to the lowest order FLR 
corrections the pressure tensor around a zeroth order double adiabatic pressure with no heat fluxes, misses the 
corrections and obtains a coefficient different both in magnitude and sign for the correction as already 

evidenced in Refs. [2l|, [2^ (see also Appendix ICl). The source of this discrepancy is in the incorrect retaining of the 
magnetosonic branch polarisation in the CGL-FLR model and in the subsequent neglect of the heat flux term, whose 
contribution needs to be retained at the same order of accuracy. 


Appendix A: Kinetic derivation of the two-flnid model with full pressure tensor dynamics 

Let us multiply the Vlasov equation for the particle distribution function /“(x, v, 7), with a = e,i denoting electrons 
and ions, by a function x(x, v). Integration over dPv, under the proper convergence conditions as |v| —>■ oo, leads to 
the so-called “phase-space conservation theorem”, 


f) 77 *^ 

+ V, • n“(v x)“ - (v • V,x)“ - — 
ot 


(F“ • V„x)“ = 0, 


(Al) 
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where F“(x,l) = ( 7 “(E(x, i)+v/cx B(x, i)) and we have introduced the particle density n°‘{x,t) = J f°‘{x,v,t)cPva.nd 
the mean value over the a-species distribution (A(x, t))“ = (l/n“(x, t)) f A(x, v, t) /“(x, v, t) By introducing 
the mean velocity u“(x, t) = (v)“, the continuity equation lEa. (IA2ll l and the Euler equation (Eq.(jA3|) are obtained 
from the zeroth (y = m°‘) and from the first order moment (y = m°‘Vi) respectively. The evolution (Eg. llAdp i of 
the pressure tensor n“ (x, t) = n°‘m°‘{viVj)°‘ — is obtained from the the second anisotropic moment, 

Xij = rn^'viVj. This equation involves the gradients of the heat flux tensor Q‘^-j = — u'l){vi — uf){vj — "«“))• 

No explicit dependence on the electric field is present since the two terms /m°‘{EiU^ + Ejuf) can be re-expressed 
in terms of the other variables by using the momentum equation. Adopting a Cartesian tensor notation with lower 
indices for the spatial components, the three moment equations for the species a are: 


dn°^ 

~dt 


dx. 




) = 0 , 


(A2) 


dt ^ dxk 


q 

~ (cE /2 Sum'll EjYi'j 

me 


1 dUf, 

dxk 


(A3) 


dt 


dQ 


kij 


dxk 




ttQ: 

k 


^ + dXk^^^ 


+ 


dxk 


n: 


ik 


(A4) 


+ SjlmBfiBm) — 0 . 


In the latter Ejim is the Levi-Civita symbol in three dimensions. Provided some closure condition for Qfjk given, 
the system of fluid equations above is closed once it is coupled to the equations for the e.m. fields, 


dxi 


— {n^q^+n^q^), 
47r 



(A5) 


dB, _ dEk 

dt ~ 


(A6) 


_ 4^ 1 dE, 

^ijk ~r 


dxj 


dt ’ 


Ji = rfq'^ul + rdq'u]. 


Notice that the set of Eas. (IA2IIA7l) leads to an energy conservation equation of the form 

tr{n“}' 


5 


E 


(m“)" + 


Stt Stt 


(A7) 


(AS) 


-V- E 


Q“ -f u“ • n“ + u° 


/tr{n“} 


^a^ocf„,ot\2 


(it“)" 


2 2 

where on the r.h.s. the heat flow vector Qf = Qf-f^Sjk/^ has been introduced. 




Appendix B: Two-fluid double adiabatic CGL closure from the full pressure tensor equation 

The CGL closure for the species a is obtained from Ea. (IA4l) . in the limit of a sufficiently strong magnetic field 
and/or sufficiently weak velocity strain, by performing an expansion in powers of 1/(rj/flo.) <C 1, where Vm“| ^ 

and fla = \q°‘\Bo/{m°'c), and by considering time scales long with respect to 11“^ such that d/dt ^ I/t)). To leading 
order we obtain 


TtO,Q: td J_ ct- R _ n 

“r c,-^m — '-'5 


0,a j 


" 3 I 


(Bl) 
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which gives 


= + (B2) 

where P|“ and P“ are the parallel and perpendicular pressures, respectively. By contracting Ea. dAdp with Sij we 
obtain 





) + (n“ % 


dxk 


2 n: 


ik 


dxk 


= 0 , 


where 


d 

dt°^ 


dt+Kdk. 


Similarly, contracting Ea. (IA4l) with bibj gives 


|(n“ 6.6,) + A «n“6.6,) - n“ I (6.6,) 


(B3) 


(B4) 


(B5) 


-\-bibj 


TT*^ 


dxk 


n 


a 

'^^dxj 


= 0 . 


Eas. (IB31IB5|) are valid regardless of the relative ordering between and and since Bm + Bm) 

disappears after contraction with either Sij or bibj. In the double-adiabatic ordering, substituting the leading order 
term for 11“ into Eas. (IB3IIB5l) allows us to derive an equation for the time evolution of P|“ and P“ . Noting 

that n^°’“^9(6i6j)/9t = 0, since |b| = 1 by definition, we obtain 


rlP° 

dt°‘ 


F)i / ^ 


(B6) 


dxk 


OXk 


(B7) 


Eas. (|B61IB7l) represent the most general writing of the double adiabatic equations, as deduced in a two fluid model. 
These reduce to the usual CGL closure [l^ of the single-fluid model when we can assume |it®| ^ |it*|. Eas. (IB61IB7l) . 
extended to include l®*-order ELR corrections and the vectorial heat-flux contribution, were recently re-derived in 
Ref. [ 2 ^, by ordering uj/^a ~ ~ Pa/Ln with Ljj^ characteristic spatial gradient of the fluid velocity. 


Appendix C: Small FLR limit for the magnetosonic waves at perpendicular propagation 


We order w/fl ~ kvth/^ ^ e <C 1. We see from Ea. dTTl) and Ea. (fT^ that the ordering Uxjxth ^ ~ e° 

and Bxy/{nomv^fj ^ e can be consistently assumed while from Ea. dTSl) we can assume Hyy/(riomv^iJ ^ e°. However, 
we see from Eg. (1121) that the ordering Uyjvth ~ £ follows, which is indeed coherent with the LE magnetosonic 
branch polarisation, for which Uy/ux ~ *£ ('Ea. (E^ '). This latter condition implies through Ea. dHl) the ordering 
{Byy — B-xx)/^ . We may now write 

Byy = PI + n w , flxx = PI- nW , (Cl) 

where the apex “(1)” labels the ^ correction. This is to ensure the conservation of the trace of Hy at any order 
in e, as suggested by the sum of Ea. dOl) and Eg. dTSl) . 

Using Eg. (IB7l) . that gives P° = 2P^kux/uj, from Eo. dT^ we obtain (cf. Eq.(A17) of [ 2 ^) 


B-xy 


i 


P° kux 

~2 IT' 


(C2) 
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Here, only the ^ e contribution has been retained, the other terms being at least of order ~ On the contrary the 
last term in Eg. (ITU) leads to an ^ contribution: combining Ea. (IC2|) with Ea. dCll) we obtain 




_n(i) 

yy 


—i 


kUy 


UJ ku^ 


(C3) 


This last equation differs from the usual CGL-FLR contribution because of the last term (see Eq.(A16) of [2^ for 
dy = 0). This term is neglected when the maximal ordering Ilyy — £ related to the assumption Uy/vth 

Ux/vth ^ is performed in Ea. flS . as is normally done in the CGL-FLR model. Solving the linear system of 
Eas. (llHll2ll41fT5]) by using Eas. (ICHIC3l) leads to the dispersion relation and group velocity of the small FLR limit 
of the LFB lEas. (13211^ 1. which gives an opposite dispersive correction with respect to the CGL-FLR magnetosonic 
mode fEas. ((Mll55)) L besides having a further contribution linear in 

Though Eas. (l32ll551) are consistent with the appropriate limit of Ea. (l551) when the heat-flux contribution is neglected, 
this latter assumption in a full kinetic treatment implies the restriction to small values of (3. In this regard we recall 
that it has been pointed out in Ref. [2l| how to recover also in the fluid /3 ~ 1 regime the correct kinetic small-FLR 
limit of the magnetosonic dispersion relation. 

At /3 ~ 1 the ~ k'^v^f ^/contributions to the small frequency dispersion relation should be retained. For 
consistency we must retain also the ~ contribution to the n = 0 harmonic in K 22 of Eq. (SOI) arising 

from the next order term in the expansion of the Bessel function, which leads to Eas. (l54ll55l) and then to the kinetic 
dispersion relation of the magnetosonic branch that includes the FLR correction to the sound term, too (EqlSni The 
kinetic dispersion relation (1551) differs from Eq. (15^ because of the factor 5 in the term. The discrepancy 

is due to the neglect of the gradient of the heat flux from Eq. ([3]) , whose lowest order contribution to the linear system 
was shown in Ref.[2l| to enter with a term dxQ xxi, ^ in Ea. (IT4l) . This thermal flux contribution, first derived 
in Ref. [ 3 ^ and more recently re-discussed in Ref. [3^, would correct the r.h.s. term of Ea. (IG3l) with a further term 
—ikQxxy/{2fl), which modifies Ea. (IG3l) into 



_n(i) 

yy 


. P^ kUy 


UJ P± kux 


l-2k^pl — 

UJ^ 


(G4) 


The above discussion is summarised by the long wave-length limit of the dispersion matrix m for the magnetosonic 
branch, as obtained by using both Eas. (IG3IIC4)) . 




[M] = 


2 I 2 

Ca + Vth 


1 - 


k^Pt\ , k^p: 


2^2-| 


ikpi 


kV: 


th 


\ 


— ikp; 


kvth 


(G5) 


In the standard GGL-FLR approximation the kpi terms are inconsistently set to zero everywhere but in the off- 
diagonal components. Including the k'^pf/8 term in square brackets of the Mxx element, which is due to the last 
r.h.s term of Ea. (IG3l) . makes it possible instead to recover the correct magnetosonic dispersion relation of Eg. (153)) in 
the small-,5 limit. If we let kdi ~ kpi, i.e. /3 1, the lowest order contribution of the heat-flux gradient should also 

be retained, and from the further contribution of Eo. (IG4l) to Eo. (IG3l) the — A:^p^X/4 correction in the Mx^x element 
follows. In this case Eg. (1551) is recovered. 

We conclude by noting that, while the assumption Uy/vth ^ £° is incorrect to describe the propagation of mag¬ 
netosonic waves in this specific geometry where Bq • fe = 0, it becomes legitimate when the condition Bq • fc = 0 is 
relaxed or when compressionless equilibrium conditions dependent e.g. simply on u = (0,My(x),0) are searched for, 
as done in Ref. [ 2 ^. In the first case, a contribution proportional to k^^c^ due to the y-component of the flJ x b/(ne) 
force of Eq.([2|) would appear at r.h.s. of Eo. dT^ . which allows us in principle to order Uy/vth ^ Ux/vth ~ at 
least as long as k^^/k± is not ordered small (i.e. of order e itself). In the second case, since Uy{x) only enters in the 
incompressible non-gyrotropic equilibria devised in Ref. [ 2 ^, no relative ordering is needed between Ux and Uy (and 
heat fluxes are null at equilibrium). We finally note that corrections to the inconsistent GGL-FLR description of the 
jrr^ag ation of fast magnetosonic waves have also been discussed in the different context of the Landau-fluid models 
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